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I. INTRODUCTION 

Historically, the notion of a physical energy contained 
in a compact 3 dimensional spacetime region (quasi-local 
energy) has proved to be rather complicated to identify 
and various attempts have been made to define a physi¬ 
cally meaningful quantity for this concept in general rel¬ 
ativity (see for a review). In the presence of spheri¬ 
cal symmetry, however, (in general, without asymptotic 
flatness) the physical mass-energy is commonly identi¬ 
fied with the Misner-Sharp-Hernandez construct the 
task becoming much more involved when the spacetime 
is not spherically symmetric. There is a large consensus 
in the relativity community that, in th^eneral case, the 
Hawking-Hayward quasi-local energy [a, l7| is an appro¬ 
priate concept to define the mass contained in a compact 
region of spacetime. 

Conformal mappings of spacetimes are widely used in 
cosmology and black hole physics {e.g., [l| and references 
therein), as well as in Penrose-Carter diagrams bring¬ 
ing spatial infinity to a finite distance in asymptotically 
flat spacetimes {e.g., i), not to mention their utility 
in alternative theories of gravity {e.g., i). The con¬ 
formal transformation properties of other quasi-local en¬ 
ergy concepts have been discussed in the literature: see 
Q for the Arnowitt-Deser-Misner mass and @ for the 
Brown-York quasi-local energy. However, in spite of the 
importance of the Hawking-Hayward quasi-local energy, 
a transformation formula for this quantity under confor¬ 
mal spacetime mappings is not available in the litera¬ 
ture except for the special case of spherical symmetry, in 
which the Hawking-Hayward mass reduces to the Misner- 
Sharp-Hernandez mass [71, llOll. for which the transforma¬ 
tion was derived recently [ll|. The present paper fills 
this gap in the literature. It is necessary to first analyze 
how the various quantities appearing in the definition of 
the Hawking-Hayward mass transform under conformal 
mappings, which is done in the next section. The follow- 
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ing section derives the desired transformation property 
(eq. (1351) , which is the main result of this paper). We then 
check that, in the special case of spherical symmetry, this 
formula reproduces the known one for the transformation 
of the Misner-Sharp-Hernandez mass. 

II. HAWKING-HAYWARD QUASI-LOCAL 
ENERGY 

The Hawking-Hayward mass Mh [1,0 is a functional 
of the spacetime metric gab and of an embedded (space¬ 
like, compact, and orientable) 2-surface S defined by 

where TZ is the induced Ricci scalar on S, 9± and are 
the expansion and shear tensors of a pair of null geodesic 
congruences (outgoing and ingoing from the surface S), 
ijj°‘ is the projection onto S of the commutator of the 
null normal vectors to S, ^ is the volume 2-form on the 
surface S, and A is the area of S. 

A. The surface S 

A co-dimension 2 surface S' in a 4-dimensional space- 
time is defined as the set of points where two independent 
functions take constant values. Let us call these functions 
and Then, 

S := {a; I (/>^^^{x) = n {x \ (t)^^\x) = . (2) 

Different 2-dimensional surfaces can be defined by choos¬ 
ing different constants 4)^'^, indeed entire 2 parameter fo¬ 
liations of the 4-dimensional manifold can be defined in 
this way for certain choices of the functions in special 

circumstances. 

Even before invoking a metric, we can associate two 
1-forms with the surface S by 


la ■■= , 

ria ■= , 


( 3 ) 

( 4 ) 
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whose restriction to the set S defines the ‘normal direc¬ 
tion’ 1-forms “ since there is no metric yet we don’t have 
a notion of orthogonality and we do not have a canonical 
way to define ‘normal vectors’ by raising the indices on 
la and Ha- These 1-forms clearly satisfy 

^ [a^b] = 0 = V [a'^b] (5) 

which are the conditions for the 1-forms I and n to be 
closed, dl = 0 = dn. This implies the metric-independent 
condition of hypersurface orthogonality on both I and n 

n A dn = 0 = I A dl. (6) 


Despite appearances, this property is not sufficient to 
define the basis vectors and would only do so in 

the case that the commutator of L and N is orthogonal to 
the surface 5, that is ha;,[L,iV]“ = 0, a condition known 
as ‘surface-forming’. In general, the partial derivative 
contains also a component tangent to S: 

= L + r where hab'f’^ = 0 (13) 

and similarly for N. Such considerations play a role in 
the 2-1-2 formulation of general relativity (see e.g. [II)- 
Note that while I (and n) is tangent to an affinely 
parametrized geodesic. 


The condition of hypersurface orthogonality for a 1- 
form is more general than the closure condition: 1-forms 
proportional to the differential of a function, fdg (where 
g and / are functions), are in general not closed but are 
hypersurface orthogonal 


(fdg) A d (fdg) = fdg A df A dg = 0 . (7) 


That is, we are free to scale the 1-forms I and n by ar¬ 
bitrary functions without spoiling their hypersurface or¬ 
thogonality. 

We note that there is also considerable freedom in the 
functions defining S since any other set of two (suf¬ 
ficiently nice) independent functions will 

define the same set S. 

Introducing now a metric gab, we impose that both I 
and n be null g°‘^lalb = 0 = g°'^nani, and from them form 
the scaled 1-forms 


la := 




na ■■= 




( 8 ) 


where l-n = g°‘^lanb. In general, null normals do not have 
a canonical scaling in contrast to the timelike and space¬ 
like cases where we can normalise to =Fl. Our particular 
choice of scaling for the null normals implies that 

tria = -1 (9) 


so that the induced metric on the surface S is written as 


Kb = gab + laUb + IbUa- (10) 

We note also that in the null case there is a restricted free¬ 
dom on the functions ; only un-mixed reparametriza- 
tions '!/;0)(^(0) leave the null condition on the normals 
intact. 

One often finds in the literature the double-null con¬ 
struction based on the vectors and iV“ defined by 

70 a 

T“:=—, N^:=^ (11) 

l-n l-n 

which have the property of Lie-dragging the 3-surfaces of 
constant 


( 12 ) 


= = Iv6(u“) = 0, (m) 

T (nor n“) is not, but is still tangent to a non-affinely 
parametrized geodesic. Defining e™ := V—Z • n we have 

fyf = e-”^Z“Va (e-”^Z^) 

= e-”"Z“Z^Vae-”" -b e-2™/“VaZ^ 

= (^-TVam^ t 

= ( —£j77i) I (15) 

where C is the Lie derivative. A similar result holds for 
ri“. 

The quantities appearing in the definition of the 
Hawking-Hayward energy are defined in terms of Lie 
derivatives in the directions I and n. The expansions 
9± and shears are defined by 


0+ 

:= I K^Cj Kb , 

(16) 

6»_ 

:= I K^C-n Kb , 

(17) 

<^lb 

■— hcd hab^-\- 7 

(18) 

^ab 

:= hlhfCn hcd - habO- , 

(19) 


while the anholonomicity LOa is defined by the commuta¬ 
tor of Z“ and n“ (or 1 and n“): 

- 77 hab {I , — 77 T ~^ab R; ^ ] , (^0) 

Z I - n Z I ■ n 

where [, ] is the Lie bracket. Note that the anholonoinic- 
ity is invariant under rescaling of the normals. These 
definitions are textbook-standard except for the anholo¬ 
nomicity uja which deserves special mention. Usuallj0 
when one is dealing with a single null congruence r;“, 
one introduces an auxiliary null vector such that 
= —1 and defines a 2-form uJab called the twist 


^ For example see Refs. M or M- 


L(</,(2)) = 1 = 
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form by cVd- This object is important because it 

vanishes if and only if is hypersurface-orthogonal (it is 
proportional to the gradient of a scalar). This object is 
not independent of the auxiliary vector |13l| but its square 
uiab is. One can show that if oJab = 0 for some choice 
of then it is zero for all choices of so the vanishing 

of ujab is a ^“-independent condition |14| and hence hy¬ 
persurface orthogonality does not depend on the choice 
of auxiliary m“. What we have here is somewhat differ¬ 
ent. We are dealing here with a single 2-surface which 
is the intersection of two null 3-surfaces. The vectors 
and n“ certainly are hypersurface-orthogonal individu¬ 
ally to these null 3-surfaces and hence both possess van¬ 
ishing ujabi but this is not what the anholonomicity uja 
is measuring (otherwise it would trivially vanish in this 
2-1-2 formulation). The anholonomicity describes how the 
two and n“ vectors ‘weave together’ (integrate to) a 
2-dimensional hypersurface, whether they are really tan¬ 
gent to a genuine 2D submanifold. Now this ‘integral’ 
sub-manifold, if it exists, is not S, but the orthogonal 
complement since the 1°' and n“ vectors are orthogo¬ 
nal to S. 


III. TRANSFORMATION OF THE 
HAWKING-HAYWARD ENERGY UNDER 
CONFORMAL MAPPINGS 

We now consider the relationship between Mh defined 
with respect to two different metrics gab and gab which 
are related by a conformal factor 


Of course, since the surface S is defined independently 
of the metric, we are considering how the energy content 
contained within a single surface changes when one scales 
the metric as above. 

In the expression ([T|) (with respect to the metric g) 
of the Hawking-Hayward mass of a 2-surface S in the 
conformally scaled world, the area of this surface is 

A = [ ji= [ d?x\/}i= [ (fxVn^= [ fj.n^. (22) 

Js Js Js Js 

Unless is constant on S, one cannot extract it from the 
integral and we have 


vCI = 


'Jsj^ 




(23) 


Only when fl is constant on S can we write \/A = fl'/A. 
This expression will be used in the integral defining the 
quasi-local energy. 

We will make use in what follows of the contracted 
Gauss equation Q 

7^ + 9+9- - i (T+ = h^^h^'^Rabcd (24) 

which provides a very useful expression of a part of the 
Hawking-Hayward energy in terms of simpler quantities. 
In the conformally rescaled world we have 


7^- 


M--2- 


= h-^V”^Rabad. 


(25) 


gab — gab- 


( 21 ) 


Using the fact that h°‘^ = and the well known transformation property of the Riemann tensor under 

conformal transformations {e.g., 0, p. 466) 


Rabc'' = Rab/ + 2df,V6] Vc In H - 2g‘^^ga[a^b]^e In H + 2V[„ In H In H 

- 2V[, In In H - 2ga[aS^^g^^Ve In H V/ In H (26) 

and, lowering one index, 

Rabcd — gdsRabc — gdsRabc (^^) 

it follows that 

1 
2 

-fV[a In H Vc In H V[„ In H g^a^d In H V-'" In H V/ In H , (28) 

where is a 2-dimensional Kronecker delta. By computing the terms appearing on the right hand side of this 
equation, one obtains 


7^-f i (T++ 2H-2 /i““h^<5f„V,]V,lnH-/i““/i%,V,,]VelnH 


InUd^jVclnU - InUgbjcVdlnU = lnHVt,lnH, 

9c[a5l^ In H Vy In H . 


(29) 

(30) 

(31) 
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Putting everything together gives 

1 


7^- 


- - ^ = n-^[n + 0+e_ - ; a+ a- 


'ab^- 


+ _ab 


'ab^- 


+ (Va In rj Vb In - Va Vf, In 0) - 2^2"^ In 12 Vc In . 


Using now the fact (which we will prove below) that 


= 12 ‘^UJabO°‘ . 


(32) 

(33) 


the integrand in the Hawking-Hayward quasi-local energy is seen to transform as 


1 


7 ^ + e+e_ - - - 2 = 12-" ( + e+e. - - aj, - 2 


1 


_+ _ab 


212 ' 


2 

lab ( 2Val2Vf,12 

VaVbl2\ 

V=12Val2' 


r V 

12 y 

122 


(34) 


The Hawking-Hayward quasi-local energy itself then transforms according to 

1 I A f r,ab/2Val2V6l2 VaV^,12^ V^12Vcl2 


Mh = 


1 / 


Stt y 1071 7s 


1 


iTZ + — 2 tiJaW“ I -|- 


47r y IGtt 7s 




122 


12 


122 


= 'U"'' + iV4F// 


^ /2Val2V6l2 Va^b^\ V°12Vcl2 


122 


12 j 


122 


(35) 


r 


One can also arrive at this final result (1351) by trans¬ 
forming directly the ingredients from which Mjj is con- show that 
structed. The 1-forms la and Ua in eqs. m and o are 
defined independently of the metric and therefore are un¬ 
changed under conformal transformations. 


— ^ai — ^a* 


(36) 


As an aside, and for ease of computation, it is easy to 
t 

1 

i _ haA^bB^^ {C^hAB) 

( 41 ) 




Hence 


la ■■ = 


la 


= 12 


-gb%hc 

la 

yJ-^-'^gb<^ll,nc 

la 


whence it follows that 


e+0. - = 12-2 ( 0 n _ - + 


^ab^- 


1 


-“t 


ab^- 


+ 


yJ-g^Hbllc 


12-® (z“Val2^) (n'’V6l22) 
Jj%(n“Val2) + ^0a (fVb^) 


:= 12 la 

and, similarly, ria = ilfia- We also have 

f=12-^r, #“=12-in“. 


(37) 


(38) 


(42) 

The Ricci scalar of the 2-surface S transforms under 
the conformal transformation as 

n = n-^n-:^h‘^^DaDbn+-^h^’>iDan){Dbn ), (43) 


Note that in the conformal world we still have the con- where Da is the covariant derivative on S which is defined 

with respect to the metric hab and whose action is related 
to the 4-dimensional Va by, using the example of a (1,1)- 
(39) tensor, 


dition 

I Tie = — I 

so that the metric is still expressed as 

hab — gab 4” ^a^b 4” ^6^a ■ 


DaX^ := hth^fhlXdXi . 


(44) 


(40) We can re-express everything in terms of the 4D covariant 
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derivative as follows 

Dan = hiv^n , (45) 

DaD^n = hihiVaDdn , 

= KhiVa (46) 

so that 

h‘^^DaDbn = (K^en) 


where dfl‘^ 2 ) = ^ i® the line element on the 

unit 2-sphere and R is the areal radius. Under a con¬ 
formal transformation with conformal factor U = U(t, R) 
(to preserve spherical symmetry), this line element be¬ 
comes ds^ = ni^ds^ and the areal radius is i? = Ui?. 
The Misner-Sharp-Hernandez mass in the rescaled world 
is simply [TTj | 

r3 

Mh = nMH - — R^v'^nWcR. (5i) 


= 


-f [fud + rfld^ 




= + (fVen^ 


h’^'^VaUd 


-f (r VeU) 


= -h {f^en}j % -h (n^VeU) Oj . 

(47) 

In the last line we have made use of the equivalent defi¬ 
nition for the expansion 0^ := h°‘^VaXi-b and similarly for 
9j. We are then left with 

= (V„U) (Vf,U) 

- ^ (r VeU) On- ^ (r VeU) 0j . (48) 

Note the cancellation of the cross terms in Oj and 6n when 
combining eqs. (l48l) and (1421) . 

We also have 


~dab 


21 ■ n 

i— [n-^ln-^nY 


2 U H ■ n 

1 


2U- 


Kb [l,n]°' = Wa 


(49) 


In spherical symmetry, a sphere S of constant radius 
R has area A = AttR"^, \f\ = f2-\/A, and h°‘^ = 
diag (i?“^, sin“^ d). In order to check whether 
eq. (I55l) reproduces eq. (I5T|) correctly, one needs to com¬ 
pute the various terms in the integrand of eq. (I55|) . The 
integrand is constant over a sphere of constant radius and 
the integral reduces to the product of the integrand and 
the area 47ri?^ of this sphere. We have 

= nUii , (52) 


and 

n f ( 2V‘=uVcU\ 

Stt V igtt Js^\ n'^ ) 


R^ 

-—v^nvan. (53) 


Since U = 0(t, R) we have VaU = U Sao + n' 6ai, where 
a dot and a prime denote differentiation with respect to 
t and R, respectively. Then it is a^lVbn = 0 and we 

are left only with the quantity —to compute. 
We have 


VaVbO = K^dadbn - K^Vabdcn 


= K^deen + K^d^^n - {K^v‘22 + Kn 


where we have used D = K^lb = U ‘^g°‘^lb- Combining 
these terms we arrive at the result (1551) obtained above. 


= -T (ro,n + nx') - (rS,^ + rj.n'). 

(54) 


IV. SPHERICAL SYMMETRY 

As a check of the formulae derived above, it is useful 
to discuss spherically symmetric situations. In spheri¬ 
cal symmetry the Hawking-Hayward quasi-local energy 
reduces [3, Il3 to the Misner-Sharp-Hernandez mass 
0. The transformation property of the Misner-Sharp- 
Hernandez mass under conformal mappings was recently 
reported in Ref. [HI- A spherically symmetric metric can 
be put in the form 


Using the line element (l50)) . the Christoffel symbols are 
easily computed: 


r^3 = o, 


(55) 

R 1 

Rsin^d 

(56) 

H ’ ^ 33 - 

B ’ 


which yields 

20' 

K^Va^bn = — (57) 
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and 


R 

Sir 



n 

R 




OTT 

R^n' 


B 


(58) 


This is the third term on the right hand side of eq. (EU; 
in fact, 


td2 

- i?2 V"OVe-R = -R'^g’^^Va^ Sbi = -— (59) 

B 

and the transformation property m of the Misner- 
Sharp-Hernandez mass is indeed a special case of eq. dSSI). 


V. DISCUSSION 

Contrary to other notions of quasi-local energy, the 
transformation property of the Hawking-Hayward quasi¬ 
local energy is not discussed in the literature. This 
concept requires one to analyze the two null (outgoing 
and ingoing) geodesic congruences associated with a 2- 
surface S, and to compute their optical scalars and an- 
holonomicity. We have presented the relevant calcula¬ 
tions, taking care to impose the normalization of the 


null normals used in the Hawking-Hayward quasi-local 
energy definition i 0, in both the original and the 
conformally rescaled worlds. Putting together the var¬ 
ious pieces, one arrives at the desired transformation 
property (l35l) for Mh- The result is reproduced by the 
calculation which involves the contracted Gauss equa¬ 
tion. In the special case of spherical symmetry, in which 
the Hawking-Hayward quasi-local energy reduces to the 
Misner-Sharp-Hernandez mass , we recover a pre¬ 

vious result on the transformation of this quantity under 
conformal rescalings m further confirming our general 
result. 

Given the wide use of conformal transformations in 
cosmology and black hole physics, we expect the result 
obtained here to be useful in these areas. An application 
to cosmological perturbations will be reported elsewhere 
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